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Abstract 

Non-adiabatic transitions are studied in a spin-boson model with multiple scattering 
points. In order to generalize the Landau-Zener formula, which describes the case 
of a single scattering point, we define an "effective gap" for a set of scattering 
points. The generalized formula agrees very well with numerical results of the non- 
adiabatic dynamics, which we obtained by a direct numerical method. This will 
make the Landau-Zener formula yet more useful in analyzing experimental data of 
magnetic-moment inversion. 
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1 Introduction 



Macroscopic quantum tunneling [1-4] has been an intriguing subject for sev- 
eral years. A few experiments have been carried out [1,4] for mesoscopic mag- 
netic systems in order to determine the coupling matrix element (or the gap) 
between the two macroscopic states. The basic assumption in analysis of the 
experimental data is the system Hamiltonian of the form 

n = \h5a^ - ^Ha^, (1.1) 

where the up-spin state | |) and the down-spin state | |) denote the two 
macroscopic states, while 6{> 0) in the first term denotes the coupling matrix 
element between the two states. The coupling hS is also called the gap, the 
reason of which is obvious in the schematic energy spectrum in Fig. 1. 

The major experimental method of measuring the gap hS is to apply an os- 
cillating field to the system and seek resonance due to the gap at if = 0. 
However, this method may pick resonance of various other degrees of freedom 
in real materials and mix them up with the resonance due to the gap [2]. 
Hence, one of the present authors [5] suggested applying a constantly chang- 
ing field instead of an oscillatory field and estimating the coupling h6 by the 
use of the Landau-Zener formula [7,8]. 

Suppose that the system is first under a large field in the negative direction 
[H <^ —h6) and the spin is approximately in the down-spin state. If the field 
is switched to the positive direction very slowly, or adiabatically, the system 
follows the ground state and eventually converges to the up-spin state for 
H ^ h5; thus the spin is flipped. If the fleld is changed more rapidly, however, 
the system does not necessarily follow the ground state and may remain in the 
down-spin branch even after the field becomes positive. In other words, the 
system is excited to the higher energy level at H = with a certain probability 
p. This phenomenon is called non-adiabatic transition. Hereafter, we refer to 
the point where the transition takes place as the scattering point. In this case, 
the scattering point is H — 0. 

The Landau-Zener theory [7,8] gives the "passing" probability p in the fol- 
lowing form: 

(1.2) 

Thus the spin is fiipped with the probability 1 — p. Here AE is the gap at the 
scattering point; hence AE — h5 in the case of Eq. (1.1). The parameter v 



p — exp 
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denotes the rate of the field change; that is, the field is increased from a large- 
magnitude negative value to a large positive value as H{t) = vt, where t is the 
time. Hence, to measure the hipping probability under a constantly changing 
field should give an estimate of the gap h5 from Eq. (1.2) [5], provided that 
the above Hamiltonian (1.1) approximates experimental systems well. 

In real materials, however, there must be always coupling between the mag- 
netic moment and phonons. Hence, in the present paper, we study non-adiabatic 
transitions in the following spin-boson model [9] : 

n = ~hSa^ - \H{t)a' + Ua + I) + EK<^'(^1 + ba). (1.3) 



The third term denotes various phonon modes and the fourth term denotes 
the coupling between the spin and the phonons. As we show below, the energy 
spectrum of the model (1.3) have many scattering points at various values of 
H. The main purpose of the present paper is to show how we should modify the 
Landau-Zener theory for such a problem. To generalize the formula (1.2) to 
multiple-scattering cases, we introduce a new parameter, the effective gap. The 
effective gap, when it replaces the gap AE in the Landau-Zener formula (1.2), 
explains well the multiple-scattering data which we numerically calculated 
using the exponential product formula [10,11]. 

If the scattering points are well separated from each other in the energy spec- 
trum, a non-adiabatic transition at each scattering point can be regarded as 
an independent event and should be well approximated by the Landau-Zener 
formula (1.2). We can explain the entire time evolution as a simple sequence of 
these independent events. This is the case if the model (1.3) has only one boson 
mode. We first consider this case in Sec. 2 and show that the Landau-Zener 
formula is indeed consistent with the numerical results. 

However, this naive description should fail if some of the scattering points 
are close to each other. This situation arises when we consider multiple boson 
modes. We show in Sec. 3 that the numerical results in two-mode cases can 
be well parametrized with a single parameter, namely, the effective gap. 

For numerical calculations of the dynamics of the system, we changed the 
magnetic field from — i/max to ifmax in the time range — tmax < ^ < ^max; hence 
the rate of the field change is = i^max/^max and the change schedule is 

H{t)^vt^^^t for - i^ax < i < imax- (1-4) 

Imax 

We made the maximum field strength i^max much greater than h6 to ensure 
that the numerical results well approximate the non-adiabatic dynamics with 
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-f^max = oo. We calculated the time evolution of the system, taking the ground 
state for H = — i/max as the initial state. Further computational details are 
given in Appendix. 

Incidentally, wc hereafter use dimensionless parameters for simplicity, taking 
-f^max as the unit of energy. A tilde mark on top of each parameter indicates 
that the parameter is dimensionless: e.g. Ti. = Ti/Z/max, ^ = ^^/-f^max, = 
huJa/H^a^, Va = Va/H^a^, H (t) = H {t) / H^^, and i = tH^^/h. Hence the 
Hamiltonian (1.3) now reads 

n = -\5a^ - ^H{i)^' + E'^a + I) +j:VaCT%bi + h^) (1.5) 

with the field-change schedule given by 

Th 1 

H{t) = Vt with V = — 2 ^ = r for - *max < * < imax- (1-6) 

max ''max 



In this dimensionless notation, the maximum field is i^max — 1- The Landau- 
Zener formula reads 



p — exp 



[1.7) 



2 Non-Adiabatic Transition and the Landau— Zener Formula 



In the present section, we treat the spin-boson model with a single boson 
mode. We calculate time evolution from the initial state and analyze it with 
the Landau-Zener formula. 

The energy spectrum of the Hamiltonian (1.5) is exemplified in Fig. 2, showing 
its dependence on H. (Note that the spectrum is calculated for each fixed 
value of H, not under a changing field H{t).) The values of the Hamiltonian 
parameters are fixed to 



S = 0.05, cDi = 0.3, and Vi = 




throughout this section. (The reason of the seemingly strange choice of Vi is 
that the spin-boson coupling often scales as l/^S^.) 

In Fig. 2, the states with positive slopes have a large component of the spin- 
down state, while those with negative slopes have a large component of the 
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spin-up state. The states with the same spin but different energies differ in the 
approximate number of excited bosons. (These interpretations are indicated 
in Fig. 2.) The initial state, which is the ground state for H = — i^max, or 
^ = — 1, is hence approximately a spin-down state without any bosons. As H 
is increased at a finite rate v, the state evolves along the ground state of the 
Hamiltonian before it encounters a nearly degenerate state at the scattering 
point at H — 0. Non-adiabatic transitions occur afterward. 

For example, the magnetization of the spin (cr^)j and the total energy E(t) — 
(7Y)f evolve as shown in Fig. 3 when the magnetic field is increased as scheduled 
in Eq. (1.6) with v = 1/500. Here the angle brackets (■ • •)f denote the expec- 
tation value with respect to the wave function at time t: (^'(t) | ■ ■ ■ | ^(i))- 
(See Appendix for details of the computation.) The initial state (approxi- 
mately I i, 0)) is scattered when it passes the scattering point at H — 0. The 
resulting state for < t < 0.3 tmax is a superposition of the ground state (ap- 
proximately I t)0)) ^iid the first-excited state (approximately | i,0)). Hence 
the magnetization is not fully inverted at H = 0. (The oscillation of the mag- 
netization after the non-adiabatic transition is due to the spin precession. This 
becomes evident when we see that the energy of the spin does not change as 
rapidly as the magnetization does; see Fig. 4.) 

The above explanation is supported by the time-evolution data of the eigen- 
state weights p,, where each pi denotes the weight of the {i — l)th excited 
state in the evolving state ^{i). In other words, we define the weights by 
Pi{t) — 1^1(^)1^, where we expand the evolving state ^(i) with respect to the 
eigenfunctions {(f)i{H(i))} of the Hamiltonian at time i as follows: 

^{i)=J2Mi)MHii)). (2.2) 



The result of the calculation for the present case is shown in Fig. 5. The state 
^{i) for < t < 0.3 tmax is the superposition of the two states | t,0) and 
I i, 0), as is expected. 

As the apphed field is increased further from H — 0, the component p2 of the 
evolving state ^(f) undergoes a subsequent non-adiabatic transition around 
t = 0.3 tinax and spawns the component ps. The component P2(t) for 0.3 t^ax ^ 
t < 0.6 tmax corresponds to | f, 1), a spin-up state with one boson emitted, 
while the spawned component pa (t) in the same time range is approximately | J, 
, 0) . The component pa undergoes further non-adiabatic transitions afterward. 
Thus the final state for f = 00 will be a superposition of spin-up states with 
various numbers of bosons. 

In this one-mode case, the scattering points are far apart from each other. 
We then expect that the whole event is regarded as a series of independent 
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non-adiabatic transitions, each of which can be treated by the Landau-Zener 
formula (1.7). Hence the passing probabihty after the transition at the scat- 
tering point H = should be 



P2{i) = exp 



1 



for 0<t<0.3t 



(2.3) 



while the probability that the state remains in the lower level is pi = 1 — p2. 
Here AEi corresponds to the energy gap between the ground and the first 
excited states at if = 0. Similarly for the subsequent non-adiabatic transition 
around i — 0.3 imax, the passing probability is given by exp[—^7i{AE2yv~^], 
where AE2 is the energy gap between the first and second excited states at 
the scattering point. The passing probability for the set of the two scattering 
points with the gaps AEi and AE2 should be given by 



P3(t) = exp 



--n{AE,rv-' 



exp 



--n{AE2fv-' 



= exp<; --TT 



{AE,f + {AE2f 



for 0.3 imax < t < 0.6 t^^. 



(2.4) 



On the other hand, the probability of the first excited component in this time 
range is 



P2(^) = exp 



■lAAE^rv-' 



1 — exp 



^7i{AE2yv-' 



for 0.3 imax < t < 0.6 t^^. 



(2.5) 



We can continue the calculation for subsequent transitions. 

Indeed, the above theoretical argument explains our numerical data of the 
magnetization very well, as is shown in Fig. 6. Thus we can in turn estimate 
the energy gap of each scattering point from the time dependence of the mag- 
netization [5,6]. 

However, we cannot naively apply the Landau-Zener theory to cases where 
some of the scattering points lie closely to each other. This is the central issue 
of Sec. 3. Before we proceed to the next section, a remark on Eq. (2.4) is in 
order; the expression can be formally regarded as if the system underwent one 
non-adiabatic transition at a single scattering point with an "effective gap" 
AE(^*^) defined in 



(2.6) 
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as has been pointed out by Kayanuma [12]. We show in the next section that 
the passing probabihty can be written in the Landau-Zener form with an 
effective gap even when scattering points are close to each other, ahhough the 
value of the effective gap can be different from the one given in Eq. (2.6). 



3 Successive Non-Adiabatic Transitions and Effective Gap 

In the present section, we consider the spin-boson model with two boson 
modes. We discuss the case where non-adiabatic transitions cannot be re- 
garded as independent events. We introduce the effective gap to generalize 
the Landau-Zener formula to such a case. We study how the magnitude of the 
effective gap depends on the distance between two scattering points. 

In this section we keep the first boson mode the same as in Eq. (2.1) and 
change the parameters of the second boson mode as follows: 

(i) 0)2 = 0.35 and % = 0.1/V(l35; 

(ii) (Da = 0.31 and V2 = 0.1/V(l3T; 

(iii) 0)2 = a)i = 0.3 and V2^Vi^ 0.1/ VO. 

We calculated the system dynamics, again taking the ground state for if = —1 
as the initial condition. 

In the case (i), there are two scattering points quite close but still apart. We 
show that the passing probability for this case can be given by the product of 
two passing probabilities with the real values of the energy gaps. 

This is not true in the case (ii), where the two scattering points are much closer 
than in the case (i). However, we show numerically that the same formula of 
the passing probability is still applicable, although we have to replace the real 
values of the gaps with an effective gap. (In fact, the real energy gaps are not 
clearly defined in the case (ii) anymore.) The remarkable point is that the 
single parameter of the effective gap determines the passing probabilities for 
various field-change rates. This is a non-trivial observation. 

In the case (iii), the two modes are degenerate and the two scattering points in 
question are reduced to one point. Then the original single-scattering Landau- 
Zener formula becomes applicable again, but with a single real energy gap 
instead of two gaps. 

We then show in the last subsection that the effective gap estimated in the 
case (ii) smoothly interpolates the real energy gaps in the cases (i) and (iii). 
Thus we conclude that the effective gap is a very useful concept in analyzing 
numerical and experimental data. 
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3.1 Case (i) 



We show in Fig. 7 the energy spectrum of the spin-boson model with two 
boson modes in the case (\): 5 — 0.05, ui — 0.3, — 0.35, V\ — 0.1/-\/(X3, 
and V2 = O.I/VO.35. The time evolution starts with the ground state at 
H = —1 and goes through the scattering points near H = 0, 0.3, 0.35,. . . . We 
calculated the time evolution of the magnetization (cr^)j of the spin and the 
total energy, as shown in Fig. 8. 

Let us focus on the scattering points near H = 0.3 and 0.35. Although these 
two are fairly separated, we cannot resolve the magnetization jump around 
H = 0.3 ~ 0.35 into two non-adiabatic transitions. The spin-precession oscil- 
lation probably smeared out the two magnetization jumps. Still, we can resolve 
it into two in the calculations of the eigenstate weights, the same quantities 
as plotted in Fig. 5. In Fig. 9, we can separate the two "scattering regions," 
where the eigenstate weights pi drastically change. Indeed, the values of the 
eigenstate weights are consistent with the Landau-Zener formulas separately 
applied to the two scattering points; i.e. the eigenstate weight p2 converges to 



P2 = exp 



1 — exp 



(3.1) 



and ps to 



P3 = exp 



X 



exp 



exp 



-^7i{AE2fv-' 



--TTiAE^fv-' 



(3.2) 



for t = tmax- Here AEi, AE2, and Ai?3 arc the energy gaps at the scattering 
points near H = 0, 0.3, and 0.35, respectively. 

Since the magnetization behaves as if there is a single non-adiabatic transition 
around H = 0.3 ~ 0.35, the situation remarked at the end of the previous 
section now becomes a reality in terms of the magnetization. That is, the 
magnetization jump around H — 0.3 ~ 0.35 is approximately described by a 
single Landau-Zener formula (1.7) with the effective gap 



(3.3) 



This is plotted in Fig. 8 as a dotted line with steps. 
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3.2 Case (ii) 



If two scattering points lie close within the width of the scattering region, we 
cannot regard two non-adiabatic transitions as independent events anymore. 
In the energy spectrum in the case (ii), the scattering points near H = 0.3 
and 0.31 are so close that we cannot even define the energy gaps AE2 and 
AE^; see Fig. 10. The time evolution of the eigenstate weights pi are shown in 
Fig. 11. We see from this figure that the two non-adiabatic processes are not 
separate enough to be independent. 

In this case, we propose a new definition of the effective gap for the non- 
adiabatic transitions around ^ = 0.3 ~ 0.31. With this newly defined effective 
gap AE'^^^\ the magnetization jump around ^ = 0.3 ~ 0.31 is described by 
a single Landau-Zener formula; see Fig. 12. 

Let us describe the definition of the effective gap in the following. First, we 
define the effective gap AE^^^^ for each scattering point with 



Pi (^max) 



i-1 



1 — exp 



--niAE, 



{eff)^2~-l 



(3.4) 



Specifically, we obtain the first estimate AE[^^^ from the value of pi (tmax) as 

AEf^ = 



2v 

l0g[l -pi(tin£«)]. 

7T 



(3.5) 



Since the scattering point at = is far apart from the other scattering 
points, the estimate AE[^^^ simply coincides with the real gap AEi between 
the ground and first excited states at ^ = 0. The second estimate AE^^^^ is 
given by 



P2(4iax) = 1 -Pl(^max) 



exp 



-IniAEf^fv 



(3.6) 



or 



AEf> 



\ 



2v 



TT 



log 



1 - 



P2(^max) 
1 -Pl(^max). 



(3.7) 



This effective gap, however, has no counterpart to be compared with, since we 
cannot properly define the energy gap between the first and second excited 
states at the scattering point near H — 0.3. This is also the case for the third 
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estimate AE^^^\ By adding the last two effective gaps squared, we tfien define 
tlie combined effective gap AE^^^^ for tfie transitions around ^ = 0.3 ~ 0.31: 

(a^(^^))' = {AEf^f + {AEf^y. (3.8) 

The dotted hne with steps in Fig. 12 was drawn, using the Landau-Zener 
formula with the above combined effective gap. 

The key point here is that thus-defined effective gaps barely depend on the 
field-change rate v; see Fig. 13. (Otherwise the definition of the effective gaps 
would be just a transformation of the eigenstate weights.) Therefore, once we 
estimate the effective gaps for one particular value of ti, we can predict the 
magnetization jump for arbitrary values of v. This is the central point of the 
present paper. 

3.3 Case (iii) 

The energy spectrum in the case (iii), uji = UJ2 = 0.3, is shown in Fig. 14. 
The two scattering points around H — 0.3 are now reduced to one. In this 
case, the second excited state at the scattering point near H — 0.3 is not 
coupled to either the first or the third excited states [12]. This is demonstrated 
in Fig. 15; the component of the second excited state, J93, entirely vanishes. 
Since the second excited state becomes irrelevant near H = 0.3, the non- 
adiabatic transition there occurs only between the first and the third excited 
states. Thus we recover the original Landau-Zener problem for two states. The 
time dependence of the magnetization is well explained by the Landau-Zener 
formula with the irrelevant states neglected. The magnetization jump around 
i7 = 0.3 is given by the formula (1.7) with the real energy gap between the 
first and the third excited states (the dotted hne with steps in Fig. 16). 

3.4 Effective gap and separation of two scattering points 

We have shown above that the case (i) and the case (iii) are explained by the 
original Landau-Zener theory for a single non-adiabatic transition between 
two states. The passing probabilities are fully described by the formula (1.7) 
with the real values of the energy gaps at the scattering points. In the in- 
termediate region including the case (ii), on the other hand, we need to use 
the effective gap for the formula, since the real energy gaps are not clearly 
defined. We here show that the effective gap estimated in the intermediate 
region smoothly interpolates the real energy gaps in the limiting cases (i) and 
(in). 
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We calculated the effective gap in the intermediate region, fixing uji — 0.3 and 
changing 002 in the range 0.3 < C02 < 0.35. (We changed V2 at the same time 
so that V2 oc 1/ \/uj2 ) By varying the value of a;2, we change the separation of 
the two scattering points around H ~ 0.3. Figure 17 shows the a)2-dependence 
of the effective gaps A^{'^\ l^Ef\ and i^Ef\ which are defined in Eq. 
(3.4). We can see that the effective gaps in the intermediate region smoothly 
interpolate the two cases (i) and (iii). This shows that the data analysis with 
the effective gap is very stable. 

To put the above in another perspective, we show in Fig. 18 the discrepancy 

between the real energy gap and the effective energy gap. In the case (i), the 
effective gap A^*^*^^) is approximately given by the real energy gaps as in Eq. 
(3.3), or 

(A^(cff) j 2 ^ 1^^^^ j 2 ^ ^^^^^ 2 ^^^^ 

We can see in Fig. 18 that the discrepancy at 0)2 = 0.35 is only about 2%. 
Although the real energy gaps may not be clearly defined in the intermediate 
region 0.3 < cj2 < 0.35, let us, only for illustration purposes, define (a) A£'2 
as the energy difference between the first and the second excited states at 
H = uji = 0.3 and (b) AE^ as the energy difference between the second and 
the third excited states a.t H = 002. The right-hand side of Eq. (3.9) thus 
estimated (the dashed fine in Fig. 18) becomes apart from the effective gap 
when the two scattering points get close to each other as 0)2 — > cui. Note that 
it is the effective gap, not the real gaps, that describes the dynamics correctly. 

In the case (iii), the real gaps tentatively defined as above is now related to 
the effective energy gap as 

A^(efr) ^ ^ /^^^ rpj^g left-hand side squared 

is the dotted line in Fig. 18. This becomes apart from the effective gap when 
the two scattering points get separated from each other. Figure 18 also shows 
that the effective gap is a useful way of interpolating the two limiting cases. 



4 Summciry 

In the present paper, we studied the time evolution of a spin-boson model 
which may be relevant to systems with macroscopic quantum tunneling. We 
introduced the definition of the effective gap to explain multiple non-adiabatic 
transitions in complicated energy spectra. With the use of the effective gap, 
the Landau-Zener formula, originally obtained for a single scattering point, 
explains the model's multiple magnetization jumps. The important point is 
that the effective gap estimated in one case of time evolution (i.e. for one 
value of v) predicts the dynamics for all other values of v. This fact will make 



11 



the Landau-Zener formula yet more useful in analyzing experimental data of 
magnetic-moment inversion. 



We remark on a related work on multiple non-adiabatic transitions. Kayanuma 
and Fukuchi [12] studied successive non-adiabatic transitions in the following 
generalized case which consists of assembly of the Landau-Zener scattering 
points: 



vt Jl J2 J3 

Jl £1 

J2 £2 

J3 £3 

; 







(4.1) 



Here the parameters v, {Si}, and {Jj} are arbitrary constants. The theory 
concludes that the overall passing probability is given by 




for arbitrary values of {Si}- 



The energy spectrum of the above Hamiltonian appears quite similar to the 
energy spectra of the present spin-boson model, particularly when ^vt is sub- 
tracted from all the diagonal elements of the Hamiltonian (4.1). In fact, con- 
sider the case where all the off-diagonal couplings { Jj} are much smaller than 
the level spacings of {si}. The structure of the energy spectrum is almost 
equivalent to the energy spectrum in Fig. 2 except for the levels | l,n > 0) 
(which do not contribute to the dynamics in any case). Then the correspon- 
dence AEi = 2Ji becomes evident. 

To put it in other words, comparison of our conclusion with Eq. (4.2) leads us 
to the speculation that our spin-boson model may be reduced to an effective 
Hamiltonian of the form 



^(eff) 



H{t) 



(eff) 



(eflf) 



£1 







£2 







(4.3) 
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However, the reduction procedure, if any, is yet to be known. 

Finally, the exponential product formula reviewed in Appendix has proven 
very useful in calculating non-adiabatic dynamics accurately. We hope the 
formula plays an important role for many problems of quantum dynamics in 
systems with more complicated interactions. 



A Exponential Product Formula for Time-Dependent Hamiltoni- 
ans 



In this Appendix we review the computational method employed in the present 
paper, namely the exponential product formula [10]. 

In our previous study [11], we discussed the time evolution of the spin-boson 
model (1.3) in a constant magnetic field using the exponential product for- 
mula. In this case, the diagonalization, if possible, of the Hamiltonian might 
suffice to know the dynamics. In non-adiabatic transitions, on the other hand, 
diagonalization of the Hamiltonian does not provide enough information for 
computation of its quantum dynamics, because of the time dependence of the 
Hamiltonian. Since quantities change dramatically when the system under- 
goes non-adiabatic transition, high precision is required to calculations of the 
time evolution; even small computational errors near the transition points can 
result in large errors at later stages of the time evolution. The exponential 
product formula is indispensable in such problems. 

When the Hamiltonian depends on time explicitly, a formal solution of the 
Schrodinger equation 

ih^nt) = nmit) (A.i) 



is given by ^'(i) = 0)^(0), where U{t,0) is the following time-ordered 
exponential: 
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When we break up the present spin-boson Hamiltonian into three parts as 



n{t) = ns{t) + nB + ni 



(A.3) 



with 



the second-order approximant U2{t + At,t) to U{t + At,t) for a small time 
step At is given [10] by 



U2{t + At,t) 

i^t^, f 1 ^ 
Hs [t + -At 



= exp 



2n 



X exp 



lAt. 



2h 



Hb exp 



exp 

iAt 



iAt 



Hb^ exp (^-^-^n 



I.e., 



U{t + At, t) = U2{t + At, t) + 0{{Atf 



(A.5) 



(A.6) 



Note that the partial Hamiltonian TYs is evaluated at time t + ^At: this is 
essential in constructing the second-order approximant to the time-ordered 
exponential [10]. Matrix elements of the approximate time-evolution operator 
are calculated in the same way as in the previous paper [11]. 
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Figure captions 



Figure 1: A schematic of the energy spectrum of the two- level system (1.1) 
(the solid line). The dashed lines denote the eigenvalues of the second term 
only. The labels | t) and | I) indicate the states in the hmit H ±00. 

Figure 2: The energy spectrum in the one-mode case (2.1). The abscissa is the 
magnetic field H. The labels | i,0) and | t^^) indicate the relevant physical 
interpretation of the states, where n is the number of excited bosons. 

Figure 3: The time evolution of the magnetization (ct^)j (dotted line) and the 
total energy E{i) (dashed line). The field was changed as H — (l/500)t. The 
ground-state energy for each value of H is also shown as the sohd line. 

Figure 4: The time evolution of the magnetization (dashed line) and the 
energy of the spin, -Espin(^) = (~^^cr^ ~ \H{t)(J^)i (solid line). 

Figure 5: The time evolution of the eigenstate weights pi. The solid line pi 
indicates the component of the ground state. 

Figure 6: The expectation value of the magnetization calculated from the 
Landau-Zener formula (solid line) for the case (2.1) with tmax = 500. The 
result is consistent with the numerical data of the magnetization if we average 
the data over the oscillations due to spin precession. 

Figure 7: The energy spectrum of the Hamiltonian with two boson modes: 
Case (i) uji — 0.3 and UJ2 — 0.35. The abscissa is the magnetic field H. 

Figure 8: The time evolution of the magnetization (cr^)f (dotted hue) and the 
total energy E(i) (dashed line) in the case (i) with v = 1/500. The ground 
state energy for each value of H is also shown as a solid line. 

Figure 9: The time evolution of the probabilities pi in the case (i) with v — 
1/500. 

Figure 10: The energy spectrum of the Hamiltonian with two boson modes: 
Case (ii) Ui = 0.3 and 0)2 = 0.31. 

Figure 11: The time evolution of the eigenstate weights Pi in the case (ii) with 
V = 1/500. 

Figure 12: The time evolution of the magnetization {a^)i (dotted hue) and the 
total energy E{i) (dashed line) in the case (ii) with v = 1/500. The ground 
state energy for each value of H is also shown as a solid line. 

Figure 13: The estimates of the effective gaps AE'^^^^ (solid hue), AE^^^^ 
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(dashed line), and AE^^^^ (dotted line) for different values of v. The estimates 
barely depend on v. 

Figure 14: The energy spectrum of the Hamiltonian with two boson modes: 
Case (iii) cui — 002 — 0.3. 

Figure 15: The time evolution of the eigenstate weights Pi in the case (iii) with 
V = 1/500. The weights pa, p^, Pq, ps, and pg vanish for all i. 

Figure 16: The time evolution of the magnetization (o"^)f (dotted line) and the 
total energy (dashed line) in the case (iii) with v = 1/500. The ground state 
energy for each value of H is also shown as a solid line. 

Figure 17: The (D2- dependence of the effective gaps. The difference between 
the solid line and the dashed line is [A^j^^''^'']^ and that between the dashed 
line and the dotted line is [AE^f^^]'^. The combined energy gap for the two 
scattering points around H ~ 0.3 is given by Eq. (3.8), or the the dashed line. 

Figure 18: The effective gap for the scattering points around H = 0.3 and 
the real energy gaps near those scattering points. The solid line denotes the 
combined effective gap (AE^'^^^)'^, the dashed line denotes the left-hand side 
of Eq. (3.3), (AEs)^ + (AEs)^, and the dotted line denotes {AE2 + AE^f. 
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Fig. 1. A schematic of the energy spectrum of the two-level system 

Fig. 2. The energy spectrum in the one-mode case (2.1). The abscissa is the magnetic 
field H. The labels | 1,0) and \ },n) indicate the relevant physical interpretation of 
the states, where n is the number of excited bosons. 

Fig. 3. The time evolution of the magnetization {(y^)i (dotted line) and the total 
energy E{t) (dashed line). The field was changed as if = (l/500)t. The ground-state 
energy for each value of H is also shown as the solid line. 

Fig. 4. The time evolution of the magnetization ((T^)^ (dashed line) and the energy 
of the spin, £'spin(i) = {—\5a^ — \H{t)a^)i (solid line). 

Fig. 5. The time evolution of the eigenstate weights pi. The solid line pi indicates 
the component of the ground state. 

Fig. 6. The expectation value of the magnetization calculated from the Lan- 

dau-Zener formula (solid line) for the case (2.1) with tmax = 500. The result is 
consistent with the numerical data of the magnetization if we average the data over 
the oscillations due to spin precession. 

Fig. 7. The energy spectrum of the Hamiltonian with two boson modes: Case (i) 
uji = 0.3 and 0)2 = 0.35. The abscissa is the magnetic field H. 

Fig. 8. The time evolution of the magnetization (cr^)j (dotted line) and the total 
energy E{t) (dashed line) in the case (i) with v = 1/500. The ground state energy 
for each value of H is also shown as a solid line. 

Fig. 9. The time evolution of the probabilities pi in the case (i) with v = 1/500. 

Fig. 10. The energy spectrum of the Hamiltonian with two boson modes: Case (ii) 
uji = 0.3 and 0)2 = 0.31. 

Fig. 11. The time evolution of the eigenstate weights pi in the case (ii) with 

V = 1/500. 

Fig. 12. The time evolution of the magnetization {(J^)i (dotted line) and the total 
energy E{t) (dashed line) in the case (ii) with v = 1/500. The ground state energy 
for each value of H is also shown as a solid line. 

Fig. 13. The estimates of the effective gaps A^(^ff), A^l^*^^ and A^J^*^^ for different 
values of v. The estimates barely depend on v. 

Fig. 14. The energy spectrum of the Hamiltonian with two boson modes: Case (iii) 
OJi = (j02 = 0.3. 

Fig. 15. The time evolution of the eigenstate weights pi in the case (iii) with 

V = 1/500. The weights ps, p^, pq, ps, and pg vanish for all t. 
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Fig. 16. The time evolution of the magnetization {(J^)i (dotted hne) and the total 
energy (dashed line) in the case (iii) with v = 1/500. The ground state energy for 
each value of H is also shown as a solid line. 



Fig. 17. The a}2-dependence of the effective gaps. 



Fig. 18. The effective gap for the scattering points around H = 0.3 and the real 
energy gaps near those scattering points. 
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